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Abstract. I review the current status of Fermi acceleration theory at relativistic 
shocks. I first discuss the relativistic shock jump conditions, then describe the non- 
relativistic Fermi mechanism and the differences introduced by relativistic flows. I 
present numerical calculations of the accelerated particle spectrum, and examine the 
maximum energy attainable by this process. I briefly consider the minimum energy for 
Fermi acceleration, and a possible electron pre-acceleration mechanism. 

1 Introduction and Motivation 

A ubiquitous feature of astrophysical objects involving relativistic flows, such 
as active galactic nuclei (AGNs), gamma-ray bursts (GRBs) and Crab- like su- 
pernova remnants (SNRs), is the presence of nonthermal, power- law emission 
spectra (i.e. with flux density F^, cx j^^", where v is the frequency and a the 
spectral index), in particular in the radio and hard X-ray or gamma-ray domains. 
This emission is believed to be produced by accelerated particles having a corre- 
sponding power-law energy spectrum; more specifically, in most of these objects 
the emission is thought to be from accelerated electrons radiating via the syn- 
chrotron or inverse Compton mechanisms (see Mastichiadis, this volume). The 
aim of the present review will be to discuss the probable mechanism of this 
acceleration and the spectra that may be expected theoretically. 

The most widely invoked mechanism for the acceleration of particles to 
power-law spectra in non- relativistic contexts, such as SNR blast waves or in- 
terplanetary shocks, is Fermi acceleration. It seems likely that shocks are re- 
sponsible for particle acceleration in relativistic flows as well, and this is indeed 
explicitly assumed in models of GRBs and Crab- like SNRs. It is then natural 
to consider how the Fermi mechanism could operate at relativistic shocks, and 
what the resulting spectrum would be. The focus of this contribution will thus 
be the relativistic version of Fermi shock acceleration. 

This review is organised as follows: in Sect. ||, I discuss the shock jump 
conditions at relativistic shocks, emphasising the aspects relevant to particle ac- 
celeration; this section is intended to be self-contained. In Sect. ^, I describe the 
Fermi acceleration mechanism in detail, first reviewing its main features in the 
context of non- relativistic shocks, and then presenting the resulting spectrum for 
ultra-relativistic and more moderately relativistic shocks. In Sect. |4| I examine 
the acceleration time scale and the maximum energy attainable by this mech- 
anism, and consider the minimum energy for Fermi acceleration of electrons in 
an electron-ion shock, and a possible pre-acceleration mechanism. 
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2 Relativistic Shocks 

The properties of shocks most important for Fermi-type particle acceleration 
are the velocities of the shock relative to the upstream and downstream frames. 
These are obtained through the shock jump conditions. 



2.1 Relativistic Shock Jump Conditions 

The shock jump conditions are derived from the laws of conservation of particle 
number, energy, and momentum. For relativistic fluids, these are, in order: 

= r2f32n2 , (1) 
r2/3i(£i+pi) =r|/32(e2+P2) , (2) 
ri(3l{ei +pi) +pi = r|/3|(£2 +P2) +P2 , (3) 

where subscripts 1 and 2 respectively refer to the upstream and downstream 
regions, n, e and p arc the fluid number density, energy density and pressure, 
all measured in the local fluid rest frame, /3 is the fluid velocity in units of the 
speed of light c, and F the corresponding Lorentz factor. The fluid velocities are 
measured in the shock frame, where the shock is stationary and both velocity 
vectors lie along the shock normal. Equivalently, (3i and f32 may be viewed as 
the shock velocity with respect to the upstream and downstream fluids. 

For simplicity, I restrict my attention in this review to unmagnetised shocks. 
The shock jump conditions for relativistic magneto-hydrodynamics (MHD) were 



reviewed by Kirk and Duffy |20 . They can yield shock jump conditions which 



differ significantly from those derived here when the magnetisation parameter, 

' - (4) 



47r(ei +pi) ' 

where B is the magnetic field measured in the local fluid rest frame, is not 
negligibly smaller than unity 



2.2 Ultra-Relativistic Shocks 

Much of the discussion of particle acceleration below will be specialised to ultra- 
relativistic shocks, i.e. those in which the shock Lorentz factor F^h = A ^ 1, and 
/3i « 1. In that case the pressure term pi may be neglected in (||) relative to the 
first term. The particles downstream of such a shock must be heated to highly 
relativistic temperatures; assuming they obey the ultra-relativistic gas equation 
of state, £2 = 3p2, one may then solve and (|^) to obtain the ultra-relativistic 
shock jump conditions, yielding the downstream velocity 



(5) 
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Also of interest is the relative velocity /3roi of the upstream and downstream 
fluids. Using the relativistic velocity addition formula, 

= ' 

the associated Lorentz factor in the ultra-relativistic limit is /^ei ~ ^sh/\/2- 

It should be noted that these shock jump conditions are independent of 
the upstream equation of state, and depend only on the downstream gas being 
ultra-relativistically hot. It will be seen below that particle acceleration in the 
ultra-relativistic shock regime thus mirrors some of the simplicity of the non- 
relativistic, strong shock regime, due to the existence of this single, well-defined 
asymptotic value (for weakly magnetised shocks) of the shock velocity ratio. 



2.3 Moderately Relativistic Shocks 

For more general values of the shock Lorentz factor .Tsh, the shock jump condi- 
tions depend on the equation of state and the temperature of the upstream gas. 
For illustration, the two opposite extremes of an ultra-relativistically hot and a 
cold gas upstream will be examined. 



Shocks in an Ultra-Relativistic Gas 

I first assume that the gas upstream of the shock already has a highly relativistic 
temperature; this might be the case for an internal shock in a GRB fireball, for 
instance, if it propagates in a medium already heated by previous shell collisions 
(see e.g. Sari & Galama, this volume). In this case, both the upstream and the 
downstream media can be assumed to follow the ultra-relativistic equation of 
state, e = ip. Equations (^ and (|^) are then readily solved to yield the jump 
condition for shocks propagating in an ultra-relativistic gas: 

/3i/?2 = \- (7) 

This relation holds for shocks of any strength, provided only that the up- 
stream gas is ultra-relativistic. The only requirement for the existence such of a 
shock solution is that the upstream flow velocity be larger than the sound speed 
in the upstream gas, which is c/ -\/3 for an ultra-relativistic gas. 



Strong Shocks and the Synge Equation of State 

I now consider the case of a strong shock, i.e. one in which the thermal energy 
upstream is negligible with respect to the bulk flow kinetic energy, so that one 
may neglect the upstream pressure pi and write ei ~ nimc^ in (|^) and (^), 
where m is the mass of individual gas particles, assumed for simplicity to belong 
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to a single species. For the downstream equation of state, I use that of an ideal 
gas of arbitrary temperature, as given by Synge ||2^: 

62 + P2 = n2mc^G i— j . (8) 

Here T2 is the downstream gas temperature, and the function G{^) is defined in 
terms of modified Bessel functions of the first kind, G'(^) = K3{^)/ K2{£,), and 
has the asymptotic expansions: 

jjic \ 41 mc _ / mc \ ^2 



Using the fact that the gas always obeys the ideal gas law p = riT, it is readily 
seen that these two asymptotes correspond to the familiar equations of state for 
non-relativistic and ultra-relativistic ideal gases, namely e = nmc^ + 3p/2 and 
e = 3p, respectively. 

With the above assumptions, the shock jump conditions (|l|-||) may be solved 
by first using (]l|) to rewrite (^) and (||) in terms of the normalised quantities 



n2mc'^ t, 

P2 _ 1 
n2mc^ ^ 



£"2 ^ = GiO - 7 ' (11) 



p, ^ = - , (12) 



where the reciprocal temperature of the downstream gas has been defined as 
^ = mc^ /T2. The two resulting equations may be solved to yield Fi and 7^2 in 
terms of ^ @: 

A = (£2+^2)^2 . (14) 

If one prefers to use Fi rather than F2 as the independent variable, the analytical 
equation (Q) , substituting the definitions (pl]-|T3[) , rnay be inverted numerically. 
The shock velocity ratio /3i//?2 resulting froni (|l3]-p^ is shown in Fig. Q as a 
function of the upstream four- velocity PiFi, along with the corresponding result 
for the ultra-relativistic upstream gas case (0). 



Electron— Ion Plasmas: The shock jump conditions derived above, based on 
the equation of state (||), are strictly speaking only valid for a gas composed of 
particles of a single mass m, as would be the case for instance in an electron- 
positron plasma. Equation ^ is readily generalised to a gas composed of species 
of different masses [p0[|2^, assuming they are in thermal equilibrium. However, 
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Fig. 1. The shock velocity ratio as a function of the upstream four- velocity, for the two 
extreme cases of a cold gas upstream obeying the Synge equation of state downstream 
{solid line), and of an ultra-relativistically hot gas upstream {dashed line) 



simulations of relativistic, perpendicular shocks in electron-ion plasmas |18| show 
that the species do not in fact achieve thermal equilibrium immediately behind 
the shock, but instead have distinct temperatures corresponding to the thermal- 
isation of their respective upstream bulk kinetic energy. The strong shock jump 
condition illustrated in Fig. || then also applies in this I now demonstrate. 

As remarked in ]20| , for a strong shock one may derive from the shock jump 
conditions (|l|-||) the relation 

£2 = ricip2c'^ , (15) 

where p2 is the total downstream rest mass density, and F^d the Lorcntz factor 
corresponding to the relative velocity (H). This equation shows that the down- 
stream energy density per unit mass is simply the upstream bulk flow energy of 
the particles as seen from the downstream frame. If this relation holds separately 
for each species, as is the case initially in the electron-ion shock simulations men- 
tioned above, one has for each species, using (pT|): 

GiO - ^ = ^rcl . (16) 

Thus although the temperatures of the species will in general be different, their 
normalised (reciprocal) temperatures ^ will be the same, and the shock jump 
conditions obtained above for a single particle mass will hold in this case also. 

Energy exchange between the electrons and the ions does take place down- 
stream of the shock in the above-mentioned simulations, but appears to result 
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in a power-law tail of the electron energy distribution rather than simple heat- 
ing, as discussed further in Sect. 4.4, One could envision shock jump conditions 
taking into account this phenomenon in the downstream equation of state; how- 
ever, an equally important component of more realistic shock jump conditions 
is the energy and momentum carried away by the strong electromagnetic pre- 
cursor emitted by the shock front . As neither of these two phenomena can 
be predicted quantitatively at present, it seems premature to attempt to obtain 
more accurate shock jump conditions for electron-ion plasmas than that shown 
in Fig. for a simple Synge equation of state. 



3 Fermi Acceleration and the Spectral Index 

In this section, I first review the basic ideas of the Fermi acceleration mechanism 
in the context of non-relativistic shocks, then discuss in some detail its applica- 
tion to the opposite extreme of ultra-relativistic shocks and the spectral index 
resulting in that case, before addressing the more involved intermediate case of 
moderately relativistic shocks. 



3.1 Non-Relativistic Shocks 

While the essential concepts of the acceleration mechanism that bears his name 
date back to Fermi [|l^ , their application to shock acceleration was first proposed 
in 1977-78 in four independent papers [|j3J|,^; of these, I will follow most 
closely below the treatment given by Bell \M. 

The acceleration scenario is illustrated schematically in Fig. |^: a high-energy 
particle, assumed for simplicity to be already relativistic, diffuses through the 
medium on either side of the shock by scattering on magnetic irregularities. These 
may be, for instance, Alfven waves self-consistently excited by the diffusing high- 
energy particles ||^ . Assuming the local Alfven velocity is much smaller than the 
shock velocity, to lowest order the magnetic scattering centres may be considered 
at rest with respect to the fluid, so that the scattering events do not change the 
particle energy in the local fluid rest frame. 

Consider, then, a particle diffusing upstream or downstream while preserving 
its energy in the corresponding rest frame. A particle initially having energy 
Ei in the upstream medium will eventually cross the shock, its velocity upon 
crossing making an angle 9^^ with the shock normal (see Fig. H). Its energy 
measured in the downstream frame, will then be given by the appropriate 
Lorentz transformation, and preserved while the particle is downstream. If it 
re-crosses the shock into the upstream medium, this time at an angle 9'^^, its 
final energy upstream, E{, will be given by the combination of the two Lorentz 
transformations : 

^ = r^ill - /3,elM^d)(l + AcIMLu) , (17) 

where /Sroi and -Trci sue the relative velocity of the upstream and downstream 
media and the corresponding Lorentz factor, and I have introduced the nota- 
tion /i for cos 9. Here and in what follows primed and unprimed quantities are 
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Fig. 2. Schematic representation of one cycle of shock acceleration: a shock propagates 
with velocity Psh into the undisturbed (upstream) medium to the right; the velocity of 
the shocked (downstream) medium relative to the upstream one is /3rei- A relativistic 
particle diffuses through the media on both sides of the shock, crossing and re-crossing 
it at incidence angles and d^u, with initial and final energies and E{ 

respectively measured in the downstream and upstream rest frames. The only 
approximation made in deriving ( p^ is that that the particle is highly relativis- 
tic, so that the rest mass contribution to its energy may be neglected. 

For non- relativistic shocks, prc\ 1, the angular distribution of these scat- 
tered particles crossing the shock may be approximated as isotropic, so that the 
flux-weighted averages of the direction angle cosines, over the relevant ranges 
—tt/2 < 9^^ < and < 9'_^^^ < 7r/2, are respectively (^^d) ~ —2/3 and 
(m'-*u) ~ 2/3. The average energy gain ( |l7| ) per shock crossing cycle then re- 
duces to 

(S) ^ ^^''-"^ " ^^^""^^ ^ + 1^'-^' ■ ^^^^ 

Fermi shock acceleration is sometimes referred to as the Fermi I mechanism 
because this energy gain is of first order in the velocity Prci- 

While a particle upstream will always eventually cross the shock, at least in 
the simple case considered here of an infinite, plane-parallel shock, once down- 
stream the particle has a certain probability of being advected away and never 
re-crossing the shock. This escape probability may be evaluated from the ratio of 
the average fluxes of particles escaping far downstream and crossing the shock. 
For a downstream flow velocity ^2, again assuming isotropy of the relativistic 
particle distribution, it is given by Pose — 4/32- 
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The combination of the energy gain factor (18) and escape probabihty Pcsc 
leads to a power-law spectrum in the accelerated particle energy, with a spectral 
index depending solely on the shock jump conditions: 

where r = /3i//32 is the shock velocity ratio. For a strong, non-relativistic shock 
in a monatomic gas, r — A, leading to a spectral index (r + 2)/(r — 1) = 2 for the 
accelerated particles. In reality, spectral indices somewhat steeper than this value 
are often observed; this difference may be due to the pressure of the accelerated 
particles modifying the shock structure (e.g. and references therein). 



3.2 Ultra-Relativistic Shocks 



I now turn my attention to ultra-relativistic shocks, i.e. those for which Fsh 3> 1 
so that the shock jump conditions derived in Sect. 2.2 apply. 



Energy Gain and Upstream Particle Dynamics 

For a downstream particle to cross the shock into the upstream medium, it must 
have 1 > fi'^^^ > /3gjj = i, so that the factor (1 + P,-eiiJ,'^^) in ( |l7| ) is always of 
order unity. If fi_d is approximately isotropically distributed, as might be the 
case for a population of relativistic particles already present in the undisturbed 
upstream medium, the factor (1 — /3reiM-d) is in general also of order unity. Thus 
in the first shock crossing cycle, a large initial boost in energy can be achieved, 
Ei/E{ r^^^ as envisioned in Q. 

For all subsequent shock crossing cycles, however, the distribution of //^d 
will be highly anisotropic; this is an essential difference between non-relativistic 
and relativistic Fermi shock acceleration [ p2|p7| ] . For an ultra-relativistic particle 
with Lorentz factor 7 ^ Fghj the kinematic condition to cross the shock into 
the upstream medium reduces to 9^^ < 1/Ah- As shown in |l^, for reahstic 
deflection processes upstream the particle cannot be deflected very far beyond 
this 'loss cone' before the shock overtakes it, so that 9^^ ^ iZ-^sh as well. In this 
case the energy gain factor reduces to 

El ^ 2 + (r.hg^d)^ ^ 1 + .2n^ 
Ei~ 2 + (r,h0^u)2 ^ 1 + M^d ' 

where the shock crossing cycle is now considered from downstream to upstream 
and back. 

The range of possible energy gain factors can be assessed by considering two 
opposite extremes for the upstream particle dynamics: deflection by a regular 
magnetic field and scattering by small-scale magnetic fiuctuations. In terms of 
the correlation length £ of the magnetic field, these two regimes respectively 
correspond to Ri^/Fsh ^ £ and Ri^/Fsh S> where i?L is the Larmor radius 
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of the particle. For regular deflection, it can be shown that for ingress angles 
< -Tsh^^u < 1, the egress angle satisfies 

l<^3h^^^d<2 ^ ^>M'-.d>-^, (21) 

while for direction-angle scattering, the direction angle at the average shock 
recrossing time satisfies 

(O^T^-^u- (22) 

In both cases, it may be seen that the typical energy gain AE' = Ef — E[ is thus 
of the order of E[ itself 



Numerical Calculation of the Angular Distribution 

As was seen in the case of non-relativistic shock acceleration, the power-law index 
of the accelerated particle distribution depends on the average energy gain per 
shock crossing and the escape probability. For relativistic shocks, both of these 
are strongly dependent on the angular distribution of particles crossing the shock, 
which as suggested in the previous section is in general highly anisotropic. Thus 
the quasi-isotropic approximations used in the non-relativistic case do not apply 
here, and the distribution of the shock crossing angles /i^u and /i^d has been 
evaluated numerically. 

For simplicity, I will focus in this section on the case where both the upstream 
and downstream particle dynamics are dominated by scattering of the particle 
momentum direction. In other words, it is assumed that magnetic fiuctuations 
dominate over the regular magnetic field in determining the particle transport, 
at least along the shock normal direction which is of interest here. This is a 
highly plausible assumption downstream, where shock-generated turbulence is 
likely to give rise to disordered magnetic fields significantly stronger than the 
shock-compressed upstream field, as is assumed in relativistic fireball models of 
gamma-ray bursts and their afterglows (see Sari & Galama, this volume) . 

One method of computing the accelerated particle distribution is through 
numerical simulations, and I illustrate below the results of such a calculation, 
after summarising the algorithm: since the nature of the particle transport up- 
stream and downstream is by assumption independent of particle energy, it is 
computationally more efficient to decouple the dynamical problem from the en- 
ergy gains. A numerical approximation to the function fdifJ-'^ui Z^'-^d)' distri- 
bution of downstream egress angles fj,'^^ for a given ingress angle, is thus first 
constructed by Monte-Carlo simulation of the downstream scattering process for 
a grid of /i'^^ values. The upstream dynamics are represented by a similarly ob- 
tained upstream egress angle distribution /u(/x^d; M^u), and both distributions, 
along with the energy gain formula (^), are subsequently used in a Monte-Carlo 
calculation of the steady-state fiux of accelerated particles crossing the shock. 



The results of such a calculation were summarised in |14 . The influence of 



the highly anisotropic injected particle distribution was seen to disappear at 
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Fig. 3. Asymptotic downstream angular distribution of the particles crossing the 
shock, showing both the flux F{^!) [dotted line) and density n(/i') [solid line) obtained 
by Monte-Carlo simulations, along with the density obtained by the eigenfunction 
method (dashed line). All distributions are normalised to unity 

a little more than a decade above the downstream injection energy, at which 
point the self-consistent angular distribution was established with a power-law 
in energy, F{E',fj,') oc F{fi')E'~P . Here F represents the steady-state flux of 
accelerated particles crossing the shock, per unit energy E' and direction angle 
cosine /i'. The asymptotic angular distribution obtained with this simulation 
method is displayed in Fig. |[ which shows both F{^') and the corresponding 
density distribution, n{^') oc F{fi')/{fi' — /Sgh)- The latter is compared with 
the distribution obtained with the very different semi-analytical eigenfunction 
method of Kirk et al. |2l| , showing excellent agreement between the two methods. 

Spectral Index and Comparison with Observations 

In the case considered above of isotropic scattering upstream and downstream 
of the shock due to a strongly turbulent magnetic field, a value of the spectral 
index p — 2.23 ± 0.01 is found by both Monte-Carlo simulations ||l|,|r^ and the 
semi-analytical eigenfunction method ||2l|] . For the opposite extreme in upstream 
dynamics of defiection by a regular magnetic field (see above), still assuming 
isotropic scattering downstream, Monte-Carlo simulations yield p = 2.30 Q. 
With particle transport including both a regular field and magnetic fluctuations, 
Bednarz and Ostrowski ||^ obtained in Monte-Carlo simulations values of p « 2.2 
in the limit Fsh 3> 1. A spectral index in the range p = 2.2-2.3 is thus a general 
feature of Fermi acceleration at (weakly magnetised) ultra-relativistic shocks, at 
least in the tcst-particle approximation used in all the above studies. 
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Spectral index values deduced from observations of astrophysical systems 
thought to involve ultra-relativistic shocks are consistent with these theoretical 
expectations. Early modelling of gamma-ray burst afterglow observations sug- 
gested p = 2.3±0.1[|||, and detailed analysis of the GRB 970508 afterglow 
spectrum yielded p — 2.2 [Q. While an equally detailed multi- wavelength spec- 
tral analysis has not been published for other afterglows, a value of p « 2.2 
seems compatible with most In Crab-like supernova remnants, the inferred 
spectral indices are similar: the best-fit model for the Crab Nebula spectrum 
corresponds to p in the range 2.2-2.3 The very good agreement between 
theory and observation is all the more remarkable given that in non-relativistic 
shocks, as mentioned above, the observed particle spectra often differ from the 
predictions of the simple test-particle theory. 



3.3 Moderately Relativistic Shocks 

For moderately relativistic shocks, the self-consistent shock crossing angle dis- 
tribution and the spectral index depend on the shock jump conditions assumed, 
as well as the shock Lorentz factor Fsh, which together determine the shock ve- 
locity ratio. Figure ^ shows, for illustration, the spectral indices obtained with 
the eigenfunction method of Kirk et al. |^ in the two extreme cases considered 
in Sect. 2.3, namely that where the upstream gas is 'cold' so that the shock is 
strong, and that where it is ultra-relativistically hot. 




Fig. 4. Spectral index of the Fermi-accelerated particle distribution as a function 
of the shock four-velocity, for the two extreme cases shown in Fig. ^ namely a cold 
upstream gas {solid line) and an ultra-relativistically hot one {dashed line) 
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It is readily seen that while for high Lorcntz factors {Fgh ~ 10) the values 
obtained rapidly converge to the ultra-relativistic case, for lower shock Lorentz 
factors the different shock jump conditions yield very different spectral indices. 
In particular, for moderately relativistic shocks in a relativistically hot gas, which 
might be relevant to internal shocks in gamma-ray bursts, the spectral indices 
obtained can be significantly steeper than in the ultra-relativistic case. For strong 
shocks, on the other hand, the spectral index goes smoothly from the ultra- 
relativistic value of p = 2.23 to the non- relativistic one of p = 2 as the shock 
Lorentz factor decreases. 

The above results apply to unmagnetised shocks and isotropic direction- 
angle diffusion both upstream and downstream. A non-negligible magnetisation 
parameter lowers the strong shock velocity ratio, leading to steeper values of 
the spectral index than in the corresponding unmagnetised case Q . Anisotropic 
diffusion in direction angle can have the opposite effect: calculations for the 
extreme case where the magnetic fluctuations are concentrated in the plane of 
the shock yield somewhat flatter spectra than the isotropic case, but by less 
than 0.1 in the spectral index p pT| . Anisotropic pitch-angle diffusion has been 
simulated by Bednarz and Ostrowski who obtained steeper spectral indices 
for weak scattering, the regime considered above corresponding to the limit of 
strong scattering. 



4 Maximum and Minimum Particle Energies 

I now turn to the question of the range in particle energies over which the 
spectrum derived above applies. I first discuss the acceleration time scale for 
the Fermi mechanism, then use it to derive the maximum energy attainable by 
this process at a relativistic blast wave, and consider an alternative scenario 
involving the initial boost which can reach higher energies. I also discuss the 
minimum energy for Fermi acceleration, and the need for a distinct electron 
pre-acceleration mechanism in electron-ion shocks. 



4.1 Acceleration Time Scale 

The acceleration time scale t^cc is defined as the time needed for the particle 



energy to increase by an amount of order itself. Since, as seen in Sect. B.2, this 
typically occurs every shock crossing cycle, tacc is roughly the cycle time, which 
is the sum of the upstream and downstream residence times t^p and tdn- In the 
case of deflection by a uniform magnetic field upstream, the former is of order 

1 _ E 
-Tsh^ci ~ qrshBij_c 

where q and lUc are the particle's charge and cyclotron frequency. The down- 
stream residence time depends on the downstream scattering process; assuming 
Bohm diffusion, it is roughly the downstream gyrotime, 



tup - 7T— ^ ^ „ , (23) 
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where primed quantities are measured in the downstream rest frame, as before. 
If the downstream magnetic field is simply the compressed value resulting from 
the (weakly magnetised) ultra-relativistic shock jump conditions, B2 ~ « 
y/SFshBuL, it can be shown that tdn ^ iup 0- 

Turbulence downstream may amplify the magnetic field B2 by a significant 
factor above the shock-compressed value, thereby reducing the downstream resi- 
dence time by the same factor; assuming the field reaches equipartition with the 
thermal pressure downstream, this factor will be of order c/wa, where wa is the 
upstream Alfven speed. On the other hand, in the case of scattering by small- 
scale magnetic fiuctuations, the upstream residence time is increased from the 
value ( p3| ) by a factor of order i?L/(^sh^), where ^, as before, is the correlation 
length of the magnetic field |l3) . Thus the value of tup given in (|2^) is a lower 
limit to tacc in all cases of interest. 



4.2 Relativistic Blast Waves and Ultra-High-Energy Cosmic Rays 

An immediate application of the above considerations is to the maximum en- 
ergy attainable by Fermi acceleration at the relativistic blast waves occurring in 
fireball models of gamma-ray bursts; I thus first review some basic properties of 
these models. After an acceleration stage, an initially radiation-dominated fire- 
ball enters a relativistic 'free expansion' phase, in which a blast wave is driven 
into the surrounding medium with the approximately constant Lorentz factor 
Fsh ~ V^T], where = E/{Mc^), M being the baryonic mass in which the 
fireball energy £ is initially deposited |^,^. This is followed by an adiabatic 
deceleration phase in which the blast wave Lorentz factor decreases with radius 
i?sh as Fsh oc ^sjf^^ [||. The transition between these two phases occurs around 
the deceleration radius i?dcc, given by 

where £1 « niuic^ is the energy density of the surrounding material. 

In the absence of energy loss processes, the maximum particle energy attain- 
able by Fermi acceleration is set by the requirement that tacc be shorter than 
the age of the system, which for a relativistic blast wave is simply Rgh/c. Using 
( p3| ) for tacc, the resulting maximum energy at a given blast wave radius Rsh is 

-Emax ~ qBiFshRsh ■ (26) 

Note that this is larger by a factor Fgh than a commonly used estimate resulting 
from a simple geometrical comparison of the particle Larmor radius with i?sh 
[ p!?! . This is due to features specific to particle acceleration at a relativistic blast 
wave, in particular the fact that an accelerated particle typically executes only 
a fraction ~ 1/ Fsy^ of a Larmor orbit upstream before recrossing the shock. 

The evolution of the product FghRsii with Rsh implies that the highest £'max 
is reached at the deceleration radius, Rsh ~ -Rdcc- It has the numerical value 

/c \i/3 

i^max « 5 X IQis ZB^e ( j eV , (27) 
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for particles of charge q = Ze, where B-g is the upstream magnetic field, £52 
the (isotropic) fireball energy, rj^ the initial Lorentz factor and uq the upstream 
density, respectively in units of microgauss, 10^^ erg, 10^ and cm~'^, these nor- 
malising values being those appropriate for a GRB fireball expanding into a 
generic interstellar medium [ p3[ . 

Equation ( p7| ) rules out the production of ultra-high-energy cosmic rays 
(UHECRs), with energies up to ~ 10^" eV, at the unmodified, external blast 
waves of relativistic fireballs. Scenarios which postulate a GRB origin for UHE- 
CRs [^sj^oUsil] must thus invoke some other site or mechanism to reach the 
required particle energies. The idea most often put forward is that UHECRs are 
accelerated at internal shocks, where substantially higher magnetic fields, close 
to equipartition, could be present both upstream and downstream and allow the 
Fermi mechanism to reach UHECR energies. However, the particle spectrum 
in this case would likely be steeper, as argued above, reducing the efficiency of 
UHECR production. Moreover, the important issue of the escape of these par- 
ticles from the interior of the fireball to the surrounding medium remains to be 
investigated, as they could suffer significant adiabatic losses due to the fireball 
expansion before escaping. 

Another possibility for acceleration at the blast wave is that the upstream 
magnetic field might be amplified by a large factor above its undisturbed value 
due to instabilities driven by the accelerated particles themselves, as was re- 
cently proposed in the context of supernova remnants by Bell and Lucek 
It is unclear, however, how far such instabilities would have time to develop in 
the relativistic context, given the comparatively short time before the modified 
upstream medium is overtaken by the shock. 



4.3 Fireballs in Pulsar Wind Bubbles 

An alternative scenario for UHECR acceleration is based on the observation that 



the initial boost examined in Sect. 3.2 can circumvent the age limit (pT]), as it 
involves only a downstream half-cycle. The maximum energy is then set only by 
the requirement that the downstream residence time tdn be less than the age of 
the system. If one assumes that the downstream magnetic field is turbulently 
amplified close to equipartition values, energies of order 10^" eV or more can be 
reached, provided that a population of relativistic particles with sufficient initial 
energy to be boosted into this range is present upstream |l3| . 

Galactic cosmic rays with appropriate energies are present in the interstellar 
medium, but constitute only a small fraction of the upstream energy density, 
so that only a correspondingly small fraction of the fireball energy could go 
to boost these to UHECR energies. However, a situation where the surrounding 
medium consists almost exclusively of relativistic particles of the required energy 
occurs naturally in the context of neutron star binary merger events: the close 
binary pulsar systems observed in our Galaxy, which are the progenitors of these 
merger events, all contain millisecond pulsars with characteristic spindown times 
of order 10^ yr, while their spiral-in times due to gravitational radiation are of 
order 3x10'' yr [E6| . These pulsars thus fill the surrounding space with relativistic 
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particles over the lifetime of the binary system, forming a large pulsar wind 
bubble in the interstellar medium. 

While the majority constituents of these pulsar wind bubbles will likely be 
electron-positron pairs, pulsar winds also seem to contain ions p^p^ . Scaling 
this ion component to millisecond pulsar parameters, one can show that it yields 
ions with energies ^ 10^"' eV, sufficient to be boosted to UHECR energies pro- 
vided Fsh ~ 10^. This process is now highly efficient: a large fraction of the 
fireball energy can go to boost these ions to UHECR energies. Moreover, for 
typical parameters the blast wave will decelerate within the pulsar wind bubble, 
resulting in a power-law spectrum of boosted ions, 

^f«^-, (2») 

with a lower bound of ^ 3 x 10^^ eV, compatible with the inferred UHECR 
source spectrum [|l3|. This scenario thus naturally provides for the acceleration 
of UHECRs into the required energy range, with the required spectrum, and 
with high efficiency. 



4.4 Minimum Energy and Electron Pre- Acceleration 

Returning now to the Fermi acceleration mechanism proper, it has as one of its 
requirements that accelerated particles see the shock as a sharp discontinuity, 
as its treatment in Sect. || makes clear. For this to be the case, the particle 
Larmor radius must be larger than the shock thickness, which is in turn roughly 
given by the downstream thermal ion Larmor radius. Ions can thus undergo 
Fermi acceleration when they have reached a few times their downstream thermal 
energy, but electrons in electron-ion shocks must first reach a minimum energy 

- rrcirriic' , (29) 

where mi is the ion mass, before participating in the Fermi mechanism and 
acquiring its characteristic spectral index. 

Unless one assumes that these objects involve solely electron-positron shocks, 
the presence of synchrotron-emitting, Fermi-accclcrated electrons in GRB after- 
glows and Crab-like supernova remnants thus requires an electron pre-acceler- 
ation mechanism. This mechanism must bring the electron energy from that 
resulting from randomisation of the bulk upstream energy, which as seen in 



Sect. 2^ is E[yi — /loiTOoC^, to E'^-^^, a factor of the mass ratio mi/rric higher. 
An acceleration process operating over precisely this energy range is the reso- 
nant ion cyclotron wave absorption mechanism discovered by Hoshino et al. ||l^ 
in numerical simulations of highly relativistic, electron-positron-ion shocks. 

This resonant ion cyclotron acceleration mechanism typically yields harder 
power-law spectra than those resulting from Fermi acceleration: the spectral 
indices p obtained from the simulations are generally less than 2, and a value 
as low as p = 1 is predicted in a quasi-linear, steady-state approximation 
The resulting picture for the accelerated electron (and positron) spectrum in 
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an ultra-rclativistic shock containing ions is thus of a relatively hard power-law 
spectrum at low energies, steepening to the p « 2.2 spectrum characteristic of 
Fermi acceleration at a break energy given by (p9[). This might explain the flat 
radio spectral indices of Crab-like supernova remnants, as well as the two breaks 
in the Crab Nebula spectrum between radio and X-ray frequencies, only one of 
which can be attributed to synchrotron cooling. 

5 Summary 

The shock velocity ratio r across a relativistic shock is in general a function 
of the assumed upstream temperature as well as the shock Lorentz factor Fgh, 
but it rapidly tends to the ultra-relativistic limit r = 3 for Fgh <i 10. The 
ultra-relativistic Fermi acceleration regime then mirrors some of the simplic- 
ity of the non-relativistic, strong shock regime, this asymptotic shock velocity 
ratio corresponding to an asymptotic power-law index of the accelerated parti- 
cle distribution. For the specific case of isotropic direction-angle scattering on 
both sides of the shock, this spectral index is p = 2.23 ± 0.01; more generally, 
a value of p in the range 2.2-2.3 is found under a variety of particle transport 
assumptions. These values are consistent with the observed spectra of sources 
thought to contain ultra-relativistic shocks, such as gamma-ray burst afterglows 
and Crab-like supernova remnants. For moderately relativistic shocks, the spec- 
tral index depends on the shock jump conditions as well as Fsh] in particular, 
shocks in a relativistic gas typically yield steeper spectral indices than the above 
ultra-relativistic values. 

The maximum energy £'max of the Fermi-accelerated particle distribution is 
determined by the acceleration time, which is in general set by the upstream resi- 
dence time. For acceleration at the unmodified, external blast wave of relativistic 
fireballs, this yields i?max ~ 10^^ eV for typical parameters of the surrounding 
interstellar medium, ruling out the production of ultra-high-energy cosmic rays 
in this context. If neutron star binary merger events give rise to relativistic blast 
waves with Fgh ~ 10"^, these can provide an alternative scenario for UHECR 
production: ions accelerated in the pulsar wind present before the merger can 
be boosted to energies ^ 10^" eV by the blast wave with high efhciency; deceler- 
ation of the blast wave in the pulsar wind bubble yields a spectral index p = 2 
and a typical lower cutoff around 3 x 10^^ eV. There is also a minimum energy 
for Fermi acceleration, set by the requirement that the shock thickness be small 
relative to the particle Larmor radius. In electron-ion shocks, this requires a dis- 
tinct pre-acceleration mechanism for the electrons, which could be the resonant 
ion cyclotron wave acceleration mechanism of Hoshino et al. p^ . 
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